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I. Introduction 


The conditions within astrophysical environments can be derived from 
observational data on atomic and molecular lines. For instance, the density 
and temperature of the gas are obtained from relative populations among energy 
levels. Information on populations comes about only when the correspondence 
between line strength and abundance is well determined. The conversion from 
line strength to abundance involves knowledge of meanlives and oscillator 
strengths. For many ultraviolet atomic transitions, unfortunately, the 
necessary data are either relatively imprecise or not available. 

Because of the need for more and better atomic oscillator strengths, our 

program was initiated. Through beam-foil spectroscopy, meanlives of 

ultraviolet atomic transitions are studied. In this technique, a nearly 

isotopically pure ion beam of the desired element is accelerated. The beam 
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passes through a thin carbon foil (2 /ig/cm ), where neutralization, ionization, 
and excitation take place. The dominant process depends on the energy of the 
beam. Upon exiting the foil, the decay of excited states is monitored via 
single-photon-counting techniques. The resulting decay curve yields a 
meanlife. The oscillator strength is easily obtained from the meanlife when no 
other decay channels are present. When other channels are present, additional 
measurements or theoretical calculations are performed in order to extract an 
oscillator strength. During the past year, three atomic systems have been 
studied experimentally and/or theoretically; they are Ar I, Cl I, and N II. 

The results for the first two are important for studies of interstellar space, 
while the work on N II bears on processes occurring in planetary atmospheres. 


2 



II. Accomplishments 


(i) Ar I 

The resonance lines of Ar I at 1048 and 1067 A provide the basis of our 
knowledge of argon in interstellar space. The resonance lines involve the 
transitions 3p^-3p^4s[3/2] ^ and 3p^-3p^4s[l/2] respectively. Until now, 
however, the available data on oscillator strengths for these transitions span 
an uncomfortably large range (Morton 1991). We have measured a value of 1.92 ± 
0.08 ns for the meanlife of the 1048 A line. This is the most precisely 
determined f-value for this line and it falls in the middle of the range of 
previously available results. Cascades into the upper state of Ar I A 1067 
prevented our experimentally determining a meanlife for this transition. 
However, using the singlet- triplet mixing-angle formulation developed by Curtis 
(1989; 1991), a collaborator in this program, we extracted a reliable lifetime 
for A 1067 as well. The value is 7.96 ± 0.33 ns. The corresponding oscillator 
strengths are f (1048) = 0.257 ± 0.013 and f (1067) = 0.064 ± 0.003. 

These oscillator strengths are within about 10% of the ones used by York 
and Kinahan (1979) and York (1983), who studied interstellar matter in the 
solar neighborhood, but the key point is that our results are more precise. 

Thus the conclusions reached by York and collaborators for a Vir and A Sco are 
strengthened. In particular, the argon abundance, like the abundances of 
nitrogen and oxygen, is approximately solar. As a result, a solar abundance 
for sulfur is also expected, but since the sulfur abundance appears larger than 
solar, contamination from S II associated with H II regions is likely affecting 
the determination of the sulfur abundance. Moreover, our results firmly 
establish the fact there is a trend of decreasing absorption line width with 



increasing atomic mass. The line widths indicate a temperature of about 10^ K 
for the neutral gas toward the two lightly reddened stars. This warm, neutral 
gas appears in clumps of approximately 1 pc in size and thus represents only a 
small fraction of the volume along the line of sight. Most of the interstellar 
material then is hot gas associated with the Local Bubble in which the Sun is 
immersed. 

A description of this work was recently submitted to the Astrophysical 
Journal as a collaborative effort with Don York of the University of Chicago. 
The title page and abstract are included in the Appendix. An abstract for a 
poster paper describing the meanlife determinations to be presented at the 
forthcoming annual meeting of the Division of Atomic, Molecular, and Optical 
Physics of the American Physical Society (DAMOP) is also in the Appendix. 

(ii) Cl I 

The 1088 and 1097 A lines of Cl I were studied because recent theoretical 
results are in conflict with oscillator strengths derived from astronomical 
data. The astronomical data (Jura and York 1978; Federman 1986) indicate that 
the 1088 A line is the stronger of the two, while the theoretical work (Ojha 
and Hibbert 1990) suggests the converse. Ab initio calculations were performed 
by us to examine possible causes for the discrepancies, but the complexity of 
this atomic system led to the conclusion that experimental data are needed to 
resolve the issue, as was also suggested by Keenan and Dufton (1990). 

We carefully examined the low energy beam-foil spectrum of Cl in the 
wavelength range 1078 - 1112 A by superposing many wavelength scans, each 
measured with a different foil. The spectral line at 1088.06 A is identified 
with the transition to ground from the level described in LS coupling as 
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3p^(^P)3d ^ 5/2 (Radziemski and Kaufman 1969). The line is a very prominent 
feature of the spectrum and is unblended. Decay curves for this transition 
were measured at beam energies of both 130 keV and 170 keV, and a meanlife of 
3.3 ± 0.3 ns resulted. 

On the other hand, the line at 1097.37 A, due to the transition between 
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ground and the level described as 3p ( P)3d (Radziemski and Kaufman 1969) 

is very weak; while clearly present, the line appears only as a shoulder on the 
stronger transition at 1098.07 A. We measured the ratio of the intensities for 
these two lines, R(1088/1097) , at a number of foil positions, corresponding to 
times after excitation of 0.4, 2.0, 4.0, 5.5, 7.0, and 9.5 ns. Extrapolation 
to t=0 gives a value for this ratio as the Cl atoms exit the foil of R = 9.3 ± 
1 . 1 . 

This ratio can be expressed by a product of the ratios of the respective 
transition probabilities and of the respective excited state populations as the 
atoms emerge from the foil: 

A < 3d ^5/2 ^ ^3/2) 

A(3d \ 2 - \ 2 ) N 0 (3d \ 2 ) 

A plausible assumption is that the foil excitation mechanism would populate 

these excited states nearly equally. Together with the equal statistical 

weights, this leads one to conclude from the large value of R that f ( 1088) » 

f ( 1097) . Should the levels be exactly equally populated initially, the 

2 

resulting meanlife of the level would be 31 ± 5 ns. That this is 

plausible is shown by considering the intensity of the A 1097 line as a 
function of time after excitation and comparing observations with an assumed 
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single exponential decay of meanlife 31 ns normalized to the data; agreement 
obtained is quite good, thus qualitatively confirming the assumption about the 
excitation process. Were the meanlife to have been only 3ns, the agreement 
obtained would be much less satisfactory. 

For single channel decay, these experimental meanlives correspond to 
oscillator strengths with values of f ( 1088) = 0.081 ± 0.007 and f(1097) = 
0.0090 ± 0.0015. For comparison, the respective oscillator strengths 
determined theoretically by Ojha and Hibbert (1990) are 0.016 ± 0.004 and 0.042 
± 0.010, while the astronomically derived values are f ( 1088) = 0.040 ± 0.010 
(Federman 1986) and f(1097) = 0.014 ± 0.005 (Jura and York 1978). The 
theoretical work clearly has the strengths reversed; our results confirm the 
trend seen astronomically, but the experimental and astronomical values for A 
1088 differ at the 3 a level. Since the astronomical result was tied to the 
somewhat uncertain oscillator strength for A 1347, we are attempting to 
remeasure the meanlife and f-value for A 1347 to see whether or not the 
differences may be due to normalization. 

(iii) N II 

For light atoms and singly charged ions, intercombination (spin-changing) 
lines are very weak, usually several orders of magnitude weaker than allowed 
transitions under LS coupling. Tripp et al. (1991) recently called attention 

to an exception in the N II spectrum, where an intercombination line at 748 A 

2 13 3 

(2p D2 - 2p3s P^) is comparable in strength to the allowed decays of Pj. 

These authors found the 748 A branching ratio to be about 30% (0.29 ± 0.15) and 

emphasized the significance of this fact for the modeling of plasmas such as 

the upper atmospheres of Earth and Titan. They also pointed out that this 
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effect was predicted by the extensive calculation of carbon sequence oscillator 
strengths by Fawcett (1987), and suggested that other predictions of that 
calculation should be tested. If such strong intercombination lines should be 
at all common in the extreme UV, the consequences for plasma modeling would be 
important indeed. 

We are investigating the following questions with this point in mind. Why 
is the 748 A line in N II so strong, and can the same mechanism be expected to 
act in other important cases? Can semi -empirical calculations making use of 
existing theoretical and experimental results improve our knowledge of the 
oscillator strengths for these N II transitions? Can new ab initio theoretical 
calculations give better results for N II, and also for C I and 0 III? 

With regard to the first question, the unusual situation in N II 2p3s is 
caused by configuration interaction (Cl) with the "plunging" configuration 

2s2p , which is higher in energy for C I but lower for large Z. In N II the 

13 3 

2p3s P and P are pushed together by Cl with the corresponding levels of 2p . 

3 1 3 3 

(The 2p P is above 2p3s while the 2p P is below.) Because they are so 

1 3 

close in energy, the 2p3s Pj and Pj mix with each other, causing both to have 

strong intercombination lines. This effect can be predicted easily: levels of 

different S but the same J and parity will mix when they are nearby in energy. 

There is a second reason why the N II 748 A line is so strong. The Cl 

with 2s2p^ enhances the 747 A *D-*P transition at the expense of the *S-^P 

3 

transition; the interference produced in the dipole integral by the 2p 
component of the 2p3s wavefunction is constructive for the former and 

1 3 

destructive for the latter. This provides a further enhancement of the D- P 
intercombination line, since it proceeds through the singlet component of the 
P wavefunction and is proportional to the same D- P dipole integral. Such 
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effects are not easy to compute quantitatively with ab initio calculations. A 
powerful method is that used by Fawcett (1987), which begins with Hartree-Fock, 
includes configuration interaction and relativity, and makes adjustments to fit 
observed energy levels. The results of such calculations are very valuable in 
calling attention to unusual effects, such as in the present case, but should 
be checked whenever possible because the semi-empirical adjustments are not 
always unique. 

In addressing the second question, for simple cases such as this one, the 

spin-orbit mixing can be determined quite well directly from the observed 

energy levels. Curtis (1989, 1991) showed the effectiveness of this approach. 

Published results of two convincing ab initio calculations agree that the 

oscillator strength in pure LS coupling is 0.200 ± 5% (Beck and Nicolaides 

1976; Luo and Pradhan 1989). From the empirically determined spin-orbit mixing 

13 13 

we get an oscillator strength for the D 2 * Pj transition of f( D 2 - Pj) = 0.023. 

Alternatively we could start from experiment instead of theory. 
Baudinet-Robinet et al. (1990) measured a precise meanlife of 0.267 ± 0.010 ns 
for the 2p3s *Pj level. A theoretical branching ratio R = A(747)/A( tot) is 
required to convert this into an f value, but since R is large, the conversion 
is not very sensitive to it. (Fawcett (1987) gives R=0.933; Luo and Pradhan 
(1989) give R=0.952.) The results are f( 1 D 2 - 1 P J ) = 0.177 and f( 1 D 2 - 3 P 1 ) = 
0.023. 

2 

We can apply the same method to the other 2p -2p3s transitions, but we 

have less extensive and precise data from which to work. We need the 

11 3 3 

oscillator strengths for pure LS coupling for the multiplets S- P and P- P. 
Taking these from Luo and Pradhan (1989), we get the results shown below. 
Because we have consistent input data from several convincing sources for the 
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747 and 748 A lines, we list their precision at 5%. This is not true for the 
other lines; we estimate their precision at 50%. We add the results of Fawcett 
(1987) for comparison in the following table. 


LINE 

WAVELENGTH 

gf (THIS WORK) 

gf (Fawcett) 

-Ip 

U 2 *1 

746.98 

0.885 

0.943 

Ip -^P 
U 2 *1 

748.36 

0.115 

0.121 

-*P 

b 0 *1 

858.35 

0.014 

0.042 

lg _3p 

860.19 

0.002 

0.005 

3 P 3 p 
*2' *1 

672.00 

0.097 

0.102 

3p _lp 
*2 *1 

670.88 

0.011 

0.011 

3 p 3 p 
M' 1 

671.63 

0.058 

0.059 

3p _lp 

M *1 

670.51 

0.007 

0.008 

3p 3p 

*0 M 

671.41 

0.078 

0.080 

3p _lp 
*0 *1 

670.29 

0.009 

0.010 


Except for the transitions involving *Sq, where the Cl cancellation effect 
mentioned above makes precision difficult, the agreement with Fawcett (1987) is 
remarkable. From these results, we predict the branching ratio for X 748 to be 
28%, in agreement with Tripp et al. (1991). A paper describing these results 
is in preparation. 

As for the third question, work is in progress. We are performing ab 
initio calculations on the relevant configurations in the carbon isoelectronic 
sequence using a variety of well-known programs, including MCHF (Fischer 1978), 
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MCDF (Dyall et al. 1989), and HER (Cowan 1981). This work will also be 
presented as a poster at the DAMOP meeting in May; the abstract appears in the 
Appendix. 


III. Summary 

Oscillator strengths for Ar I, Cl I, and N II have been obtained. Our 
results for Ar I AA 1048, 1067, which are the most accurate to date, are 
f ( 1048) = 0.257 ± 0.013 and f(1067) = 0.064 ± 0.003. The improved accuracy 
strengthens the conclusions of York and collaborators regarding the properties 
of the interstellar medium toward A Sco and a Vir. We have completed the first 
reported laboratory measurements of the oscillator strengths for the Cl I lines 
at 1088 and 1097 A; our results are f(1088) = 0.081 ± 0.007 and f(1097) = 
0.0090 ± 0.0015. The fact that A 1088 is the stronger line confirms the 
prediction derived from astronomical data, contrary to recent theoretical 
results. In the spectrum of N II, we have used a semiempirical treatment of 
spin-orbit interaction, in combination with ab initio theory, to derive 
oscillator strengths for UV transitions. Our results confirm and explain the 
unusual strength of the intercombination line at 748 A, for which we report gf 
= 0.115 ± 0.006. With this understanding, the plasma processes associated with 
planetary atmospheres can be studied with more confidence. 
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Appendix 


Accurate Oscillator Strengths for Ultraviolet Lines of Ar I: 
Implications for Interstellar Gas Temperatures 

S.R. Federman, D.J. Beideck, R.M. Schectman 
Department of Physics and Astronomy, University of Toledo, Toledo, OH 43606 

and 

D.G. York 

Astronomy and Astrophysics Center, University of Chicago, Chicago, IL 60637 

Abstract 

Analysis of absorption from interstellar Ar I in lightly reddened lines of 
sight provides information on the warm and hot phases of the interstellar 
medium. The details of the analysis are limited by the quality of the atomic 
data. Accurate oscillator strengths for the Ar I lines at 1048 and 1067 A and 
the astrophysical implications are presented. From lifetimes measured with 
beam-foil spectroscopy, an f-value for A1048 of 0.257 ± 0.013 is obtained. 
Through the use of a semi-empirical formalism for treating singlet- triplet 
mixing, an oscillator strength of 0.064 ± 0.003 is derived for A1067. Because 
of the accuracy for these new results, the conclusions of York and colleagues 
from spectra taken with the Copernicus satellite are strengthened. In 
particular, for interstellar gas in the solar neighborhood, argon has a solar 
abundance and the warm, neutral gas is not pervasive. 
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Abstract Submitted 
for the 1992 Annual Meeting 
of the Division of 
Atomic, Molecular and 
Optical Physics 

Suggested Poster PACS Number: 32.70 Fw 

Session: 

ATOMIC & MOLECULAR 
STRUCTURE 


Redetermination of the Meanlives of the 3p 5 4s J=1 
Levels of Neutral Argon for the Investigation of Interstellar 
Matter . R.M. SCHECTMAN, S.R. FEDERMAN and D.J. 
BEIDECK, Univ. of Toledo *. -To study the structure of nearby 
interstellar matter, we have re-examined the meanlives of the 
3p 5 4s[3/2] 1 and [1/2^ levels of neutral argon to determine f- 
values for the ground state transitions at 1048.2 and 1066. 7A. 
We measured the meanlife of the [3/2] transition by beam foil 
techniques. The value 1.92 ± 0.08 ns was obtained. For the 
[1/2] transition a value was derived from the [3/2] 
measurement by using experimental term values for the four 
levels of the 3p 5 4s configuration to compute the degree of 
mixing of the J P and 3 P levels; from the mixing, the ratio of 
the meanlives of the [1/2] and [3/2] levels can be computed. 
The value obtained for [1/2]! was 7.96 ± 0.33 ns. 
Corresponding f-values were 0.257±.013 and 0.064±.003. Our 
results agree well with the weighted mean of previous 
measurements and with theory, but are of increased precision. 
They will be used to reinterpret Copernicus Satellite Data for 
neutral argon in the solar neighborhood. 

Work supported in part by NASA .grant NAGW-2457. 
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Abstract Submitted 
for the AMO Division Meeting of the 
American Physical Society 
20-22 May 1992 

Suggested session: PACS number 31 

Atomic Spectroscopy (UV) 


Calculation of UV Oscillator Strengths in N II Using MCHF 

Wavefunctions. D.G. ELLIS, The University of Toledo. — Oscillator 

2 

strengths have been calculated for the ultraviolet transition array 2p - 2p3s in 

12 3 

the spectrum N II . Recent astrophysical, experimental, and theoretical 

works have called attention to two interesting features of these lines. First, the 

13 1 

intercombination line, D 0 - P., is exceptionally strong; second, the P 

1 1 ^ 1 l 1 

branching ratio A( D- P)/A( S- P), is unusually large. Both these effects 

reflect the strong intention of 2s^2p3s with 2s2p^; the *P terms of these 

■i 

configurations cross between N II and 0 ID, while the P cross between C I and 
N D. In the present work, two different sets of Multiconfiguration Hartree- 
Fock^ solutions have been obtained for these states. (One solution set employs 
multiple orbitals of the same type and the other does not.) These wavefunctions 
have been used to compute oscillator strengths for all the lines in the 2p - 2p3s 
transition array; the results are compared with published theoretical and 
experimental values. 

•Work supported by NASA grant NAGW-2457. 

1 T.M.Tripp et al„ ApJ. 368 , 641 (1991) 

2 Y.Baudinet-Robinet et al. Phys.Rev. 42 , 1080 (1990) 

3 B.C.Fawcett, At.DaLNuc.Dal.Tables 37 , 41 1 (1987) 

4 CJLFischer, Comp.Phys.Comm. 1 , 151 (1969); 14, 145 (1978) 
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